Supersymmetry and Unconventional Quantum Hall EfFect in Graphene 
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We present a unified description of the quantum Hall effect in graphene on the basis of the 
8-component Dirac Hamiltonian and the supersymmetric (SUSY) quantum mechanics. It is re- 
markable that the zero-energy state emerges because the Zeeman splitting is exactly as large as 
the Landau level separation, as implies that the SUSY is a good symmetry. For nonzero energy 
states, the up-spin state and the down-spin state form a supermultiplet possessing the spin SU(2) 
symmetry. We extend the Dirac Hamiltonian to include two indices j| and j[ , characterized by the 
dispersion relation E {p) ocp-'T+^i and the Berry phase 7r(j| ~ ji)- The quantized Hall conductivity 
is shown to be a^y = ±(2n-|-j| + ji)2e^ /h. 



The quantum Hall effect (QHE) is one of the most 
remarkable phenomena in condensed matter discovered 
in the last centuryP, Electrons, undergoing cy- 

clotron motion in magnetic field, fill Landau levels suc- 
cessively. Each filled energy level contributes one conduc- 
tance quantum e^/h to the Hall conductivity dxy. The 
Hall plateau develops at Uxy = v{e'^ /h), where y is the 
filling factor. It tells us how many energy levels are filled 
up. Hall plateaux have been observed at v = 1,2,3,... 
in the conventional semiconductor QHE. 

Recent experimental developments have revealed un- 
conventional QHE in graphene|J, 3, 0, 0|- The filling 
factors[3,|^,,£] form a series [Fig^b)], 




J/ = ±2, ±6, ±10, ■ 



(1) 



where the basic height in the Hall conductance step is 
4e^//i except for the first step which is just one half. In 
the bilayer graphene QHE, the series reads 



±4, ±8, ±12,. 



(2) 



where all steps have the same height Ae^/h. A most re- 
cent experiment has shown a fine structure at v = 
0, ±1, ±4 in monolayer graphene. 

In this paper, we present a unified description of the 
QHE in graphene on the basis of the 8-component Dirac 
Hamiltonian and the supersymmetric (SUSY) quantum 
mechanics Our first observation is that the Zeeman 
splitting is as large as the Landau-level separation, as 
leads to the emergence of the zero energy state. Fur- 
thermore, the spin SU(2) symmetry is exact within each 
Landau level for nonzero energy states. Our second ob- 
servation is that the above Landau-level structure is a 
manifestation of SUSY at the K and K' points. 

We first formulate the QHE in monolayer graphene. 
Then, in order to make the underlying mathematical 
structure clear, we generalize our analysis to the system 
where the dispersion relation has the form E (p) cc p^^ ^-'i 
and the Berry phase is given by 7r(j| — Jx) in the absence 
of the magnetic field. The monolayer graphene corre- 
sponds to — 1 and ji = 0, while the bilayer graphene 
to ji = 2 and ji — 0. We show that the zero-energy 



FIG. 1: (a) A schematic diagram of the low-energy disper- 
sion relation near the Dirac points (K and K' points) in the 
graphene Brillouin zone. Only two Dirac cones are inequiv- 
alent to each other, producing a 2-fold valley degeneracy in 
the band structure, (b) The quantized Hall conductivity as 
a function of the filling factor i/ in monolayer graphene. (c) 
The Landau level spectrum of monolayer graphene for elec- 
trons (red) and holes (blue) at the K and K' points. The 
arrow indicates the spin. This spectrum is a manifestation of 
SUSY. 



states have the jf-fold (j|-fold) degeneracy for up-spin 
(down-spin) electrons at the K point. As a result the 
quantized values of the Hall conductivity become 



= ± 



it + h 



4e^ 
IT' 



0,1,2, 



(3) 



We also discuss the relation of our SUSY formalism to 
the Atiyah-Singer index theorem. 

We start with the SUSY description of the monolayer 
graphene. The low-energy band structure of graphene 
can be approximated as cones located at two inequiva- 
lent Brillouin zone corners called the K and K' points 
[FigGJa)]. In these cones, the two-dimensional energy 
dispersion relation is hnear and the electron dynamics 
can be treated as 'relativistic' Dirac electrons|l^, in 
which the Fermi velocity vp of the graphene substitutes 
for the speed of light. 

Corresponding to the K point (-I-) and the K' point 



2 



(— ), we have two Dirac Hamiltonians 

= v^ia^P^ ±ayPy)+ (3mvl, (4) 
where Pi = —ihdi + eAi is the covariant momentum, and 



a, 

(7, 



I2 

-I2 



(5) 



with CT* the PauH matrix for the spin degree of freedom 
and I2 the 2x2 unit matrix. We have introduced the 
vector potential Ai , in terms of which the magnetic field 
is -B = V X A. We assume a homogeneous magnetic field 
B = (0, 0, —B) with B > 0. In the zero field case, the 
Hamiltonian Q has a linear dispersion E(p) = zLvpp for 
both spin states and for both K and K' points. The mass 
of quasiparticle excitations is zero, to = 0, in the naive 
band-structure calculation for noninteracting quasipar- 
ticles on the hexagonal lattice of graphene. However, 
it is natural to expect a nonzero excitation gap when 
the Coulomb interaction is taken into account. Hence 
we have included the effective mass term into the Dirac 
Hamiltonian I^J- Our analysis is vaHd whether m = or 
TO ^ 0. 

The Hamiltonian ^ is expressed asfl^ 



Q± 



with 



Q± = VF {cTxPx ± (^yPy ) ■ 

It is diagonalized, 



(6) 



(7) 



diag. J Q±Q± + w?v' 



2, ,4 _ 
F' 



1±Q± + m^v^ J , 
(8) 



where the negative component describes holes. 

It is convenient to distinguish electrons on the K 
and K' points by assigning the pseudospin to them. 
Namely, we call the electron on the K (K') point the 
up-pseudospin (down-pseudospin) electron. We can com- 
bining them into the 8-component Dirac Hamiltonian, 

Hj^ = vyUxP-x + VFTzayPy + (3mvl = diag.{H^ , H^) , 

(9) 

where Tz is the Pauli matrix acting on the pseudospin 
space. The pseudospin is reversed by the operation 

TxHj,Tx^diag.iH^,H+). (10) 

Furthermore, the Hamiltonian has a symmetry, 

ri/±r = -H^, (11) 

with r =diag.(CT^, cr^). It impHes that a solution j^f) to 
the Dirac equation with energy E has a particle-hole con- 
jugate partner rl^*) with energy —E. Therefore, it is 



enough to explore the energy spectrum of electrons for 
> at the K point. Nevertheless, to make the un- 
derlying physical and mathematical structure clear, we 
present all the spectrum in what follows. 
We consider the quantitv(l3| 

= 2Q±Q± = 2vl {-ihV + eAf T ^ehvla^B. (12) 

Since this has the same form as the Pauli Hamiltonian 
with the mass to* = l/4w| except for the dimension, 
we call it the Pauli Hamiltonian for brevity. The energy 
spectrum £„ of the Dirac Hamiltonian is constructed once 
we know the one En of the Pauli Hamiltonian II12II . It is 
to be noticed in II12|I that the direction of the magnetic 
field is effectively opposite at the K and K' points. 
We introduce a pair of operators 



V2h 



(Px+iPy), 



V2h 



{Px - iPy) (13) 



with the magnetic length = yh/eB. The commuta- 
tion relation [a, a^] = 1 follows from [P^j-Pj/] = 
The operator Q± is expressed as 



Al 
A 



Q- 



A 
At 



(14) 



with A — hujca, where ujc — \/2vy/Ib = v^y/^eB /h. 

Let us first make a generic argument that is valid 
for an arbitrary operator A. We focus on the Hamilto- 
nian =diag.(i?"'"t^, H'^^) related to the K point, with 
H+t ^ A'^A and = AA^ . This is a simplest ex- 

ample of the SUSY quantum mechanics where the 
superalgebra reads 



(15) 



with Q+ the supercharge. The two Hamiltonians H^^ 
and _ff+t ai-g superpartners. We consider separately the 
eigenvalue problems for the up-spin and down-spin com- 
ponents. 



(16) 



where E^l[ > E+'^^ > E^^K We assume SjJ"^ = 0, as 
implies that the SUSY is a good symmetry [li|- (This is 
indeed the case in our system with A = hjOca.) Using the 
relations 



AH+^ = AA^A = H+^A, 



(17) 



we obtain 



= AH+^\^+^) = (18a) 
H+Ut|^+i) = A^H+^\^j+^) = E+^A^\^j+^). (18b) 

If £;+T ^ 0, A\ij+^) is an eigenstate of H+K Similarly, 
if 0, At|^+^) is an eigenstate of H+\ Thus, 
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there is one-to-one correspondence between the up-spin 
eigenstate and the down-spin eigenstate for nonzero en- 
ergy states [Fig[2]. They are said to make a supermulti- 
plet, since the correspondence is made by the supercharge 
Q+. We identify the up-spin (down-spin) sector as 
the bosonic (fermionic) sector at the K point. 

It is necessary to examine the cases E'^'^ ^ and 
0, separately. \iEQ^ ^ 0, the lowest energy eigen- 

n+l) 



state of is oc A\il)^^). Then, the states |?/', 

and \'^n'^) make a supermultiplet having the same energy 
E^^ . Namely, there hold the relations [Figl^a)] 



1 



+T 



1 



(19a) 

for n > 0. On the other hand, if Eq^ = 0, for 
which A^^It/ij]"^) = additionally, there hold the relations 
[FigEtb)] 



1 



+T 



+T\ 



1 



+1 



+l\ 



(19b) 

for 71 > 1. The similar analysis is appHcable also 
to the Hamiltonian related to the K' point, where 
H- ^diag.{H-\H-^) with H'^ = A'' A and R-^ = 
AA^ . It is to be noticed that the bosonic (fermionic) 
sector is identified with the down-spin (up-spin) states 
at the K' point. The energy spectrum is 4-fold degen- 
erated except for the zero-energy state. We cannot say 
anything about the degeneracy of the zero-energy state 
by this general argument. 




"fermiOTk;" (bji "bosonic" 



■4- 



A "fermionic" 



♦ — + 



FIG. 2: The energy spectrum of the superpartner Hamiltoni- 
ans H^^ and H^^ for electrons at the K point. Two states 
on the same horizontal line have the same energy, making a 
supermultiplet, except for the zero-energy state {Eq^ = 0). 
The zero-energy state is j|-fold ( j|-fold) degenerated for up- 
spin (down-spin) electrons, (a) The case with ji = 0. We 
have j| = 1 in monolayer graphite, while ji = 2 in bilayer 
graphite, (b) The case with ji 7^ 0. 



We apply the above analysis to the monolayer graphene 
with the "Dirac" Hamiltonian ijHJ with Ijl4|l . where A = 
hvcCL and E^'^ = 0. Using the commutation relation 
[a, =1, it is trivial to see that the energy eigenvalues 



of the Dirac Hamiltonian ff^ are 



and 



7 + T 



C-+T _ c+l 



±mv-n 



(20) 



= ±nwc\jn + 1 + (77124/^2^2) (21) 

for n > 0. If 771 = 0, there exists one zero-energy state for 
up-spin electrons but not for down-spin electrons at the 
K point [FigI2Ia)]. The existence of the zero energy state 
is an intriguing property of the SUSY theory, where the 
bosonic and fermionic zero-point energies are canceled 
out ^.12]- The physical reason is that the Zeeman split- 
ting is exactly as large as the Landau level separation. 
This is a well known property^] of the Dirac electron 
in magnetic field though it is overlooked in all previous 
literature on the QHE in graphene. On the other hand, 
there exists one zero-energy state for down-spin electrons 
but not for up-spin electrons at the K' point. Recall that 
the direction of the magnetic field is effectively opposite 
at the K and K' points. 

There exists the 4-fold degeneracy in the zero-energy 
state due to electrons and holes. No plateau is made 
ai v = because the Landau level made of the zero 
energy states is half-filled. However, the degeneracy is 
removed between electrons and holes once the mass term 
is present. However small the mass tti may be, the Hall 
plateau emerges at = since holes are filled before 
electrons. Thus, by assuming m 7^ 0, we can explain a 
recent experimental data^3| at J/ = 0. 

We next investigate a general case, where A is given 

by 



At 



(22) 



with j| and ji being integers such that j| > Ji- The 
monolayer graphene is given by j-\ = 1, ji = 0, and the 
bilayer graphene by = 2, = 0. The physical meaning 
of the indices is discussed after the energy spectrum is 
constructed. 

We study the energy spectrum of the Hamiltonian 
H+ = {Q+,Q+} with ((HI and lO. Although the 
Hamiltonian is written in terms of A and , the ba- 
sic physical variable is the momentum, or equivalently, a 
and at. Thus, we consider the state 



1 



111) = ^(aT)"|0). 

VTl! 



(23) 



We make an analysis at the K point. The zero-energy 
up-spin states are given by the condition A\n) = 0. They 
are |0), |1), • • • , |j| — 1), which are degenerate in IV'cj't). 
On the other hand the zero-energy down-spin states are 
determined by requiring A'^\n) = Q. If j| = 0, there are 
no zero-energy states, and the supermultiplet is given 
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by Cnil [Figlla)]. If ix ^ 0, they are |0), |1), 
\2l — 1), which are degenerate in iV'o'^)- Consequently, 
there exists the {jy + jj^)-fold degeneracy in the zero- 
energy state for electrons at K point, corresponding to j| 
up-spin electrons and Jx down-spin electrons [FiglJIb)]. 
The similar analysis is made also at the K' point. 

We apply the above analysis to the "Dirac" Hamilto- 
nian ijHJ. The eigenvalues are easily calculated by eval- 
uating A \n] and \n). For the bosonic sector, we 
derive 



and 



for n > 1, where 



(24) 
(25) 

(26) 



For the fermionic sector, in the case of = we find 

Et^^E-^ ^e+^in+n) (27) 
for n > [FigEJa)]; in the case of ji 7^ we find 

E^^ = fo"^ = ±mvl (28) 

and 

(29) 

for n > 1 [FigEJb)], where 



e+^{n) = ±n^,J^—n-JjL + __|:. (30) 



{{n-JiVr 



In the case of m = 0, the zero energy states are 4(j| 
fold degenerated and all other states are 4-fold degener- 
ated. This energy spectrum implies that the Hall con- 
ductivity is quantized as in 

The bilayer graphene corresponds to = 2, — 0. 
The energy spectrum is obtained as 



E^^ =EJ = ±mvl 



(31) 



and 



5+1 =E-^ =±h\ujc\^n{n+l) + m^vl/h'^ojl, (32a) 
for n > 1 and 



£+J- = £-T = ±h \uj^\ + {n + I) + m?vllH^UJ^^ 

(32b) 

for 71 > 0, as agrees with the previous result We 
expect the emergence of a Hall plateau at = due to 
the electron-hole splitting also in the bilayer graphene. 



We now discuss the physical meaning of the indices j| 
and ii in Ij22p by setting m = 0. We derive the energy 
dispersion relation in the zero field (5 = 0). Setting Pi = 
Pi and parametrizing p = {p cos <j>, p sin cj)) , we obtain 



A = E{p) e^^J'T-Ji)-^, A'' = E (p) e-'^J'T-J'i)-^, (33) 



with 



(34) 



The eigenvalues of the "Dirac" Hamiltonian ^ are 
zLE (p) with the eigenstates 



1=^) ^ ( ie'O'T-Ji)'^ ) 



(35) 



at the K point. Hence, the energy dispersion relation 
is given by E {p) . The Berry phase is defined as a loop 
integration 



Fb = ^ ^ Ckdxk 



of the connection field 



It becomes 



C, = {n\^—\n) 



Fb = 7r(jt - Ji) 



(36) 



(37) 



(38) 



for the eigenstates |±). Thus the indices and Ji are 
fixed by the dispersion relation and the Berry phase in 
the zero field. 

We comment on the Witten index and the Atiyah- 
Singer index theorem. We focus on the K point. The 
Witten indexfl^ is given by 

Aw = dim [ker H+^] - dim [ker = j| - Jx, (39) 

where j| = dim [kcri?+1^] and ji = dim [keriZ+J-] are 
the numbers of zero-energy states of H~^^ and H^^, 
respectively. We have explicitly shown in the present 
model that the Witten index defined in the nonzero field 
{B ^ 0) is equal to the Berry phase defined in the zero 
field (B = 0), Fb = ttAw- This is a general property, 
as we briefly sketch. The Witten index is equal to the 
Fredholm indexfl^ 

Ap = dim [ker^] - dim [ker A^j ^ (49) 



and 



ker = ker A^^ A = ker A, (41) 



ker = ker AA'' = ker A^ . (42) 



5 



Now, according to the Atiyah-Singer index theorem, Ap 
is invariant as B ^ 0, where Ap becomes the chiral 
anomaly [l^. The chiral anomaly is given by the Berry 
phase. Consequently it follows that Fb = ttAw- 

The Witten index, j| — ji , and the degeneracy of the 
zero-energy state, j^+Ji, are different objects. Neverthe- 
less, they are identical in graphene since ji — 0- Thus, by 
combining the contributions from the K and K' points, 
the degeneracy 4j| = 4 in monolayer graphene can be 
related to the Berry phase, as pointed out in Refs.m3|. 

We have presented a unified description of the QHE in 
graphene based on the SUSY quantum mechanics. The 
key observation is that the Zeeman splitting is exactly as 
larage as the Landau level separation at the K and K' 
points in graphene. It is remarkable that its consequence 
is the exact spin SU(2) symmetry in each Landau level, 
however large the magnetic field is. Indeed, an arbitral 
linear combination of the up-spin and down-spin states 
belongs to the same Landau level. Furthermore, the elec- 
tron has the pseudospin in addition to the spin. It leads 
to the SU (4) symmetry together with SU(4) skyrmion 
excitationsyjl, as in the bilayer system of the conven- 
tional semiconductor qhe[3. 

Interactions between electrons are not taken into ac- 
count in this paper, except for a possible mass term 
which would arise from Coulomb interactions. It re- 
moves the 4-fold degeneracy in the zero-energy state of 
monolayer graphene, as explains the emergence of a Hall 
plateau at = in recent experimental dataj^. How- 
ever, it does not remove the 4-fold degeneracy in higher 
Landau levels. With respect to this degeneracy, there 
are senariosjl^, ^3 that the degeneracy is removed by 
Coulomb interaction. In particular, it is pointed out|l3 
that the pseudospin SU(2) symmetry is broken explicitity 
down to U(l)xZ2. If this is the case, the total symmetry 



must be [SU(2)]spinX[U(l)xZ2]pseudospin, since the spin 
SU(2) symmetry is exact as we have emphasized. We 
wish to make a detailed analysis of a fine structure of 
Hall plateaux based on our SUSY formaHsm in a subse- 
quent work. 



The author is grateful to Professors T. Ando, Y. lye 
and Z.F. Ezawa for fruitful discussions on the subject. 
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